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ABSTRACT 

In  this  paper,  various  contact  and  crack  problems  for  an  orthotropic  substrate, 
stiffened  by  elastic  films,  are  considered.  The  film  is  modeled  as  a  membrane  and  the 
substrate  as  an  orthotropic  half-plane  with  the  principal  axes  of  orthotropy  parallel  and 
perpendicular  to  the  boundary.  The  problem  is  formulated  in  terms  of  a  system  of 
singular  integral  equations.  Various  asymptotic  analyses  are  carried  out  in  order  to 
determine  the  nature  of  the  stress  singularities.  The  special  cases  for  which  the  solution 
is  obtained  and  results  are  provided  include  the  stiffened  half-plane  without  the  crack, 
the  cracked  half-plane  without  any  stiffeners,  the  orthotropic  half-plane  with  a  single 
stiffener  and  a  crack  emanating  from  the  endpoint  of  the  stiffener,  and  a  broken 
stiffener  on  a  cracked  half-plane.  In  each  case  the  influence  of  the  relative 
crack/stiffener  dimensions,  the  film/substrate  stiffness  ratios,  and  the  material 
orthotropy  on  the  stress  intensity  factors  is  studied  and  some  sample  results  axe  given. 

1.  Introduction 

The  problem  of  stress  concentrations  caused  by  thin,  elastic  film  overlays  on 
substrates  has  attracted  a  great  deal  of  attention  in  recent  years  because  of  its 
importance  in  a  number  of  microelectronics,  optical  and  structural  applications.  The 
edges  of  these  films  are  known  to  be  regions  of  stress  concentration  [1,  2],  the  stresses 
being  produced  by  processing,  thermal  mismatch  or  mechanical  and  thermal  loading. 
The  investigation  cf  the  related  stress  analysis  problem  is  essential  for  understanding 
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the  associated  physical  effects  such  as  dislocation,  generation,  film  debonding  and  film 
and  substrate  distortion  and  cracking.  It  is  further  known  that  especially  for 
applications  in  microelectronics  and  composites,  both  the  film  and  the  substrate 
demonstrate  significant  anisotropy  in  their  elastic  properties  [3].  This  anisotropy  must 
be  taken  into  account  while  studying  the  mechanics  of  the  film-substrate  composite 
medium. 

Among  the  early  studies  of  the  problem  of  a  film  bonded  to  an  elastic  substrate 
one  may  mention  those  by  Arutiunian[4,  5]  who  formulated  the  problem  in  terms  of  a 
singular  integral  equation.  In  [4]  the  problem  of  a  single  film  and  in  [5]  the  interaction 
of  periodically  spaced  films  was  considered.  A  simpler  formulation  of  the  same  problem 
including  the  interaction  of  multiple  stiffeners  was  considered  by  Erdogan[6].  In  the 
solutions  given  in  [4]- [6],  the  film  was  assumed  to  be  very  thin  and  hence  the  effect  of 
bending  stiffness  and  the  variation  of  stresses  through  film  thickness  was  neglected. 
Thus,  the  film  was  modeled  as  a  membrane  and  the  only  contact  stress  was  the 
interfacial  shear  stress.  In  [6]  a  simple  and  efficient  numerical  technique  incorporating 
the  correct  form  of  the  stress  singularities  was  used  to  solve  the  resulting  integral 
equation.  Once  the  shear  stress  on  the  interface  was  determined,  the  stresses  in  both 
the  film  and  the  substrate  were  easily  calculated.  It  was  shown  that  the  contact  shear 
stress  had  a  square  root  singularity  having  a  Mode-II  stress  intensity  factor.  The 
magnitude  of  this  stress  intensity  factor  was  shown  to  depend  on  a  stiffness  parameter 
that  is  a  measure  of  the  relative  stiffnesses  of  the  half-plane  and  the  film.  Erdogan  and 
Gupta(7]  considered  the  single  stiffener  problem  and  studied  the  regularity  of  the 
solution  in  terms  of  the  stiffness  parameter.  Delale  and  Erdogan[8]  used  these  solutions 
to  study  the  associated  failure  problem.  The  membrane  stiffener  problem  has  also  been 
studied  recently  by  Erdogan  and  Joseph[9]  who  gave  a  detailed  asymptotic  analysis  of 
the  singular  behavior  of  the  solution  for  an  arbitrary  membrane-substrate  contact  angle 
and  showed  that  the  power  of  the  stress  singularity  varies  between  0  and  1/2  as  the 
contact  angle  varies  between  0  and  90  degrees. 

In  the  present  study,  the  contact  stress  problem  for  elastic  films  directly  bonded 
to  an  orthotropic  substrate  containing  a  crack  is  considered.  In  general,  depending  on 
their  relative  dimensions,  the  film  and  the  substrate  can  be  modeled  as  membranes, 
plates  or  half-planes  [6,  10,  11].  The  main  thrust  of  this  paper  is  towards  applications 
where  the  dimensions  of  the  substrate  axe  considerably  greater  than  that  of  the  film  and 
hence  the  substrate  may  be  modeled  as  a  half-plane  and  the  film  as  a  membrane.  Tk^ 
effect  of  material  anisotropy  is  taken  into  account  by  assuming  both  the  film  and  the 
substrate  as  being  orthotropic.  The  general  problem  is  formulated  in  terms  of  a  system 
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of  singular  integral  equations  which  are  solved  for  various  film  and  crack  geometries.  A 
particular  emphasis  is  placed  on  the  surface  crack  in  the  substrate  initiating  at  the  film 
end. 

2.  The  General  Formulation  of  the  Problem 

Consider  the  plane  elasticity  problem  for  an. orthotropic  half-plane  stiffened  by  a 
number  of  membranes  and  containing  a  crack  (see,  for  example  Figures  1  and  2).  The 
problem  may  easily  be  formulated  by  using  the  solutions  for  a  pair  of  dislocations  in  the 
half-plane  and  a  pair  of  concentrated  forces  acting  on  the  boundary  as  the  Green’s 
functions.  Let  the  half-plane  contain  a  pair  of  dislocations  at  a  point  (zj,  t/j)  having  the 
components  of  the  Burgers  vector  fx  and  f2  parallel  and  perpendicular  to  the  boundary, 
respectively.  The  stresses  at  a  point  (z,  y)  in  the  half-plane  due  to  these  dislocations 
may  be  expressed  as 

°ii  (*,  ¥•  *»,  Vi)  =  Kn{x,  y,  Zj,  i/O  fx  +  Kn  (z,  y,  x1,  yx)  f2  , 

<*x2  ( x>  y;  *i>  Vi)  =  K2i(x>  y,  *i.  Vi)  A  +  K22  (x>  v>  *i>  v\)  A ,  (!) 

<y-i2  ( x>  y;  *1,  Vi)  =  ^31(x,  y>  *i>  Vi)  A  +  K32  (x>  v>  *1,  yJ  A » 

where  the  functions  Kt]  ( i  =  1,  2,  3,  j  =  1,2)  are  given  in  Appendix  A  (see  [12]). 

Similarly,  the  stresses  in  an  orthotropic  half-plane  due  to  the  boundary  tractions 

0i2(O,  y)  =  /3*(y-y«),  *u( 0,  y)  =  A6(y-y0)>  (2) 

are  given  by 

<7n  (*»  y;  Vo)  =  Ki3(x>  y>  y«)  A  +  Ku  (x>  y>  y»)  A. 

crl2  (x,  y,  y0)  =  K23{x,  y,  y0)  f3  +  K2A  (z,  y,  y0)  /4.  (3) 

<r22  (*,  y,  y0)  =  #»(*,  y ,  ys)  /a  +  ^  (*»  y»  y»)  A- 

where  the  functions  Ktp  (i  =  1,  2,  3  and  j  =  3,  4)  are  also  given  in  Appendix  A[13]. 

The  functions  KX]  contain  the  parameter  Ch  defined  by  (A12)  which  characterizes 
the  anisotropy  of  the  substrate.  It  was  shown  by  Chou[14]  that  Ch  has  a  range  between 
-4  and  oc.  Physically,  Ch  =  0  corresponds  to  ah  isotropic  material.  It  was  further 


3 


demonstrated  that  the  sign  of  C \  has  an  important  influence  on  the  analysis  of  stress 
fields.  It  turns  out  that  the  roots  of  the  characteristic  equation  resulting  from  the 
solution  of  the  equilibrium  equations  for  the  plane  orthotropic  medium  are  real  for  Ch  > 

0  and  complex  for  Ch  <  0.  The  formulation  considered  here  is  valid  for  all  values  of  C*. 
It  should  be  noted  that  the  constant  C\  is  essentially  the  same  as  the  orthotropy 
parameter  k  defined  in  [15]  and  [16]  ( C 4  =  2  («  — 1),  see  [17]  for  discussion).  In  this 
paper  numerical  results  are  presented  for  six  different  materials  having  the  elastic 
constants  listed  in  Table  1.  Materials  1  and  2  are  fiber-reinforced  composites.  Material 
3  is  isotropic.  Material  4  with  CK  =  -  0.0022  is  nearly  isotropic  and  is  included  largely 
to  verify  the  numerical  procedure  developed  for  orthotropic  materials.  Materials  5  and 
6  are  Silicon  and  Gallium  Arsenide,  respectively. 

Since  the  bending  stiffness  of  the  film  is  assumed  to  be  negligible,  it  may  be 
modelled  by  a  membrane  in  which  the  stress  component  crn  vanishes,  a22  is 
independent  of  the  thickness  coordinate  and  the  distributed  interfacial  shear  stress  cr12 
=  f3(y)  may  be  treated  as  a  body  force.  Thus  we  may  assume  that  in  (2)  and  (3)  f4  =  0 
and  the  crack-contact  problem  described,  for  example,  in  Fig.l  may  be  formulated  by 
considering  fx,  f2  and  /3  as  continuously  distributed  functions  that  satisfy  the 
appropriate  continuity  and  equilibrium  conditions.  The  continuity  of  displacements 
along  the  film-substrate  interface  requires  that 

dy  Vl ^  =  ~§y  y)’  a'<y<  b"  1  =  1’"’  n  (4) 

where  vY  and  v2  are  the  v-components  of  the  displacements  in  the  film  and  the  substrate 
respectively,  a ,  and  6,  refer  to  the  film  ends,  and  without  any  loss  in  generality  the 
continuity  condition  is  expressed  in  differentiated  form. 

If  we  now  designate  the  stress  state  in  the  half-plane  by  a„,  cr^,  avv  and  define 
the  coordinate  system  (s,  n)  shown  in  Fig.l,  the  remaining  conditions  refer  to  the  fact 
that  the  crack  surfaces  are  traction-free,  that  is 

(Tnn(s)  =  a  sm29  4-  <7yy  cos2#  -  criy  sin26  =  0 ,  (A  <  s  <  B),  (5) 

=  (°’yy~C7zr)  C0S $  ° zy  C0326  =  0,  (A  <  S  <  B ),  (6) 

Bv  observing  that  dv/dy  =  and  that  /3  is  the  only  applied  force  acting  on  the  film, 
from  the  equilibrium  of  the  film  it  may  easily  be  shown  that 
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(7) 


O  01  } 

^  «i(y)  = - j  /3(y0)  dy0  ,  (a,  <  y  <  &,),  i  =  1,  n, 

a. 

where  the  film  is  assumed  to  be  orthotropic  having  the  stiffness  coefficients  C*;  and  0* 
is  defined  by1 


(8) 


Similarly,  on  the  boundary  of  the  half-plane  from  (1),  (3)  and  <7^(0,  y)  =  0,  criy(0,  y)  = 
f3(y)  it  may  be  shown  that 


JL 

dy 


{a 

<*o  +  { ^31(0,  y\  *i,  yi)  h(t)  dt 

B  A  n  b,  v 

+  I  #32(0,  y;  xj,  Vx)  f2(t)  dt  +  £  \  #33(0,  y,  y0)  f3i(y0)  dy0  l, 
a  1=1  i  J 


where 


cr0  =  an(x,  ±00), 

x  —  s  cosO,  y  =  s  sinO,  x1  =  t  cos9,  y1  =  t  siw9, 
3  =  _ ^LL _ 

12  {<?,.  Cn  -(<?„)■}’ 


(9) 

(10) 
(11) 
(12) 


C,j  being  the  stiffness  coefficients  of  the  half-plane  (see  Appendix  A).  In  solving  the 
problem  it  is  assumed  that  the  substrate  is  subjected  to  uniform  tension  crJV  =  o0  away 
from  the  film-crack  region.  Other  simple  forms  of  practical  loading  may  be  uniform 
temperature  changes  and  residual  stresses  which  may  be  interpreted  in  terms  of  a 
constant  strain  mismatch  along  the  film-substrate  interface.  Observing  that  0 0  is  the 
only  external  load  amplitude,  it  may  easily  be  considered  as  being  also  a  measure  of  the 
strain  mismatch  on  the  interface^] . 

Substituting  from  (1),  (3),  (7)  and  (9)  into  (4)-(6)  we  obtain  a  system  of  integral 
equations  for  the  unknown  functions  fu  f2  and  /3,,  where  f3i  is  the  interface  shear  for  the 
tth  film,  (i  =  1,...,  n).  The  conditions  of  equilibrium  of  the  films  and  the  single¬ 
valuedness  of  the  displacements  require  that  the  solution  of  the  integral  equations  be 

and  are  the  coefficients  of  the  compliance  matrix  for  the  substrate  and  the 
film,  respectively. 
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subject  to 

\ 

I  My  a)  dVo  =  0,  t  =  1,  2,  n, 


(13) 


(14) 


3.  The  Integral  Equations 

In  materials  stiffened  by  ’  astic  membranes  the  endpoints  of  the  stiffener  axe 
known  to  bf*  points  of  singularity.  Thus,  aside  from  the  cracking  of  the  membrane  due 
to  high  tensile  stresses,  the  most  likely  modes  of  failure  in  the  bonded  medium  appear 
to  be  debonding  along  the  interface  and  crack  initiation  and  propagation  in  the 
substrate  around  the  singular  points.  With  the  membrane  model,  debonding  may  be 
interpreted  simply  as  the  shortening  of  the  film  and,  therefore,  does  not  require  any  new 
concepts  and  analysis.  Cracking  of  the  substrate  normally  requires  the  solution  of  the 
problem  for  a  crack  growing  radially  from  the  singular  point  (Fig.  1 ,  A  =  0,  8  =  8cr). 
The  direction  of  the  crack  growth  must  be  determined  by  maximizing  the  ratio  of  the 
strain  energy  release  rate  to  the  fracture  toughness.  Since  the  material  is  orthotropic,  in 
the  problem  under  consideration,  both  of  these  quantities  are  dependent  on  8  and  8cr 
and  may  not  be  determined  without  the  full  fracture  toughness  characterization  of  the 
orthotropic  medium.  However,  in  orthotropic  materials  the  principal  planes  of 
orthotropy  are  also  known  to  be  the  weak  cleavage  planes.  Furthermore,  under  the 
Mode-II  conditions  prevailing  around  the  film  end  for  —  30°  <  8  <  30°  (where  8cr  falls) 
the  8  variation  of  the  strain  energy  release  rate  is  relatively  weak.  It  is,  therefore,  safe 
to  assume  that  in  the  film  problem  considered  8  =  0  is  by  far  the  most  likely  fracture 
plane  (Fig.l). 


By  substituting  0  =  0,  from  (5),  (6)  and  (4)  it  may  then  be  shown  that 


f  m 

t 


ZTjr  dt  +  }  Nu{z,  t)  /i(*)  dt  +  |  A'late  t)  /a(t)  A  +  £  j  ^13(1,  Vo)  htiv,)  dVo 

c  c  ,=  1  a- 

=  0i(x),  (c  <  x  <  d),  (15) 


m 


-3~x  dt  +  I  N2l(x,  t)  fx{t)  dt  +  N22{x,  t )  f2{t)  dt  +  £  |  N23(x,  y0)  f3t{y0)  dy0 

c  c  ,=  *  a- 

=  g2(x),  (c  <  x  <  d),  (16) 
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i  J 

N31(y,  t)fx{t)  dt  +  |  N32(y,  t)f2{t)  dt  +  j  N^y,  y0)f3](y0)  dy0 
i  a, 

=  53(y)>  (a;  <  y  <  6;)>  (j  =  1v,  »)  (17) 


where  the  Cauchy  kernels  and  the  Fredholm  kernels  Nt]  follow  from  K,-  (t,  j  —  1,  2,  3) 
given  in  Appendix  A.  Thus  referring  to  Appendix  A  it  may  be  shown  that 


„  ,  ,  2  f(C,  +  4)  (2x+((C4+2))  \ 

ll(  ’  '  *Ck\2(z+t)  [j  +  (CK  +  2)  xt+ t2  })  ' 

(18) 

Ar12(i,  t )  =  0, 

(19) 

...  ,  2j(C4+4)vJA3 

N,3(l'!/“)  *  K.  {(£  +  V  **)’  +  CK  o?AV„}  ’ 

(20) 

N 

r*. 

II 

o 

(21) 

N22(x,  t )  —  Nxx(x,  t), 

(22) 

..  ,  ,  2j(Cl+4)IJ20A3 

^.{(j^  +  v^  +  a^vd’ 

(23) 

.  if.(C,  +  4)5  iAV 

3'{v'  ]  *{(e  +  \W)2  +  chSM}’ 

(24) 

..  ,  K,(Ct  +  4)5  A  f  y 

33  >)  *  {(<*  +  AV)’  +  Ca  ’ 

(25) 

...  ,  0'„\n(y-y.) 

N^y-y-]  ~0n  h\(C„+  4)  ’ 

(26) 

11 

sT5  q 

O 

(27) 

II 

O 

(28) 

9^y)  ,j(  CA+  4) 

(29) 

where  /f(y  —  y0)  is  the  Heaviside  function  and  the  constants  CA,  Kn, 
defined  in  Appendix  A. 

\ ,  and  Ke  axe 

( 


4.  The  Special  Cases 


In  the  special  case  of  an  orthotropic  half-plane  without  any  cracks  and  stiffened 
by  a  single  membrane  along  x  =  0,  —  a  <  y  <  a,  the  integral  equations  of  the  problem 
are  reduced  to 


1 

x 


f/3  d»<  hTlic^T)  jA  M  iy° 


A  ( 7 


o  _ 


>j(Q+4) 


=  <r0,  -  a  <  y  <a,  (30) 


subject  to 

MVo)  dy0  =  0. 


(31) 


We  now  define 

y0  =  at,y=  ar,  ^  {f}  =  11  ?  XX  W  (32) 

By  using  the  following  results  [9] 


1  Tn(t)dt 


r 

|  u<t)  rB(«) 

•i 


0,  n  =  0,  —  1  <  r  <  1, 

Un-i(r)>  n  =  1,2,...,  -l<r<l, 

~  7=7=  [r  “  ( I r  I  /r)  ^r2  ~  1  ]  .  n  =  0,l, 

Nr  - 1 

dt  =  ^i(r)  'll  -  r2,  n  >  0, 


(33) 


•i  I  r  |  >1 


(34) 


and  a  method  of  collocation,  (30)  may  be  reduced  to 

00  f  A  \!  1  -  r2 1 

E4,C/-.(r,)(l+  V  =  1.  >=1.2 . 


(35) 


where  Tn(t)  and  Un(t )  are  the  Chebychev  polynomials  of  the  first  and  second  kind, 
respectively.  The  linear  system  (35)  is  solved  by  truncating  the  series  at  n  =  N  and  by 
monitoring  the  convergence  of  the  results.  The  selection  of  the  collocation  points  r,  is, 
in  general,  arbitrary.  It  is  known  that  good  convergence  is  obtained  by  selecting  more 
point;  near  the  ends [9],  For  the  present  study,  the  collocation  points  are  selected  as 


W  =  0,  r,  =  COS0„  6 ,  = 


(2  i  —  1)  x 

TN  ’ 


(i  = 


(36) 
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From  (30)  it  may  be  seen  that  f3(y0)  has  a  form  similar  to  the  stress  field  due  to  a  crack 
under  pure  Mode- II  loading  and,  therefore,  it  is  possible  to  define  a  stress  intensity 
factor  at  the  film  edges  as  follows: 


fc;(a)  =  lim  ^2(a  -  y)  <7^(0,  y). 


(37) 


The  value  of  this  stress  intensity  factor  is  dependent  on  the  stiffness  parameter  A,  which 
is  proportional  to  the  ratio  of  the  stiffness  of  the  substrate  to  the  stiffness  of  the  film  in 
the  y  direction.  Figure  3  shows  the  variation  of  the  normalized  stress  intensity  factor  at 
the  film  edge  as  a  function  of  A,.  For  A,  =  0,  that  is,  for  an  inextensible  film,  (30)  has  a 
closed-form  solution  [6]  and  the  value  of  the  normalized  stress  intensity  factor  is  unity 
which  is  the  upper  bound.  Similarly,  from  (30)  and  (32)  it  may  be  seen  that  as  the 
stiffness  of  the  film  decreases  (i.e.,  as  A,— ►  oc),  the  shear  stress  f3  and,  consequently  the 
stress  intensity  factor  k?  tends  to  zero. 


Once  the  shear  stress  in  the  contact  region  is  determined,  the  tensile  stress  in  the 
film  may  be  obtained  from  the  following  relation 


v 


1 

h 


h(y9)  dVo 


-a 


cr'n  a 


( y/*Y  g  A^jy/a) 

»  k 


(38) 


A  typical  plot  of  this  film  stress  is  shown  in  Fig. 4.  It  is  seen  that  the  maximum  normal 
film  stress  is  at  the  center  of  the  film  i.e.  at  y  =  0.  This  implies  that  the  y  =  0  is  the 
most  likely  site  for  potential  film  failure.  Once  the  film  fails  completely  in  the  center,  a 
crack  can  propagate  either  along  the  interface  or  into  the  substrate.  The  latter  case  is  a 
strong  possibility  as  demonstrated  in  [6].  For  the  case  of  two  stiffeners  placed 
symmetrically  with  respect  to  the  x-axis,  it  was  shown  that  there  was  a  “strong” 
singularity  in  the  state  of  stress  in  the  half-plane  for  the  limiting  case  as  the  distance 
between  the  stiffeners  approaches  zero.  Recognizing  that  complete  film  rupture 
represents  this  limiting  case,  one  would  expect  that  the  crack  would  almost  certainly 
propagate  into  the  substrate  or  along  the  interfaces.  This  problem  for  the  cracked 
substrate  is  considered  in  Section  6  of  this  paper. 

The  second  special  case  which  is  of  considerable  practical  interest  is  the  case  of  a 
crack  perpendicular  to  the  boundary  of  a  symmetrically  loaded  orthotropic  half-plane 
without  any  stiffeners.  In  this  case  fx  is  the  only  unknown  function  and  the  integral 
equations  (15)-(17)  are  reduced  to 
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hit) 

t  —  x 


dt  -f- 


N{x,  t)  fx{t)  dt 


1*2 

Kn 


(To  ,  (c  <  x  <  d ), 


(39) 
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where 


Mix  t)  -  2  f  ( Q  +  4)  _  {2  x  -f  t  {CK  +  2)}  ] 

*  Ck  |2  (i  +  t)  1 12  +  (CK  +  2)  X  t  +  t2}] 

the  dislocation  density  f^t)  is  defined  by 

hit)  =  -  ft{v2(t,+  0)  -  v,(t,  -0)}, 


(40) 


(41) 


and  the  constants  Kn  and  CK  are  given  in  Appendix  A. 

Note  that  for  isotropic  materials  CK  =  0  and  from  (40)  it  may  easily  be  shown 

that 


lira  N(  x,  £,  CA) 

Q-*  o 


I  j  1  ,  2  t(x  -  £)) 

*  \(z  +  t)  +  (i  +  tf  j  ’ 


(42) 


which  is  the  known  result  for  an  isotropic  half-plane  with  a  crack  perpendicular  to  the 
boundary. 


If  c  >  0,  the  function  /t(£)  must  satisfy  the  following  single- valuedness  condition: 


The  integral  equation  (39)  may  then  be  solved  by  following  the  procedure  described 
previously  in  this  section.  After  determining  fx[t)  the  stress  intensity  factors  may  be 
defined  by  and  evaluated  from 


=  -M*™  pf  -J2(<*  “  x)  /i(*)>  (44) 

==_iiI3  pf  “  c)  hix)-  (45) 


From  (39),  (40)  and  (42),  it  may  be  observed  that  the  kernel  is  independent  of  the 
elastic  constants  in  the  isotropic  case,  and  depends  only  on  a  single  dimensionless 
material  constant  Ch  for  orthotropic  materials.  The  stiffness  constants  KJ 2  in  (39)  and 
2^/(1  -f-  k)  in  the  isotropic  case  simply  serve  as  scaling  parameters  for  the  crack  opening 
displacement.  Thus  the  stress  intensity  factors  would  be  independent  of  the  material 
constants  in  isotropic  materials  and  would  be  a  function  of  C j,  only  in  orthotropic  solids. 
The  results  for  the  internal  crack  problem  are  given  in  Table  2.  From  this  table  it  may 
be  seen  that  the  stress  intensity  factors  are  monotonically  decreasing  functions  of  Ch 
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and  the  relative  distance  {  ( d  +  c)/(d  —  c)}  from  the  surface.  As  {  (d  -f  c)/(d  —  c)}  — ► 
1,  i.e.  as  c  — ►  0,  fcj(c)  becomes  unbounded  (because  of  diminishing  ligament  size  c)  and 
kx(d)  approaches  the  value  given  b  '  the  edge  crack  solution. 

A  second  observation  that  is  important  in  this  problem  concerns  the  effect  of  a 
90-degree  rotation  in  the  axes  of  orthotropy.  From  the  definition  given  by  (A  12),  CK  is 
seen  to  be  invariant  with  respect  to  a  90-degree  material  rotation.  Hence  the  kernel  of 
the  integral  equation  and  consequently  the  stress  intensity  factors  would  also  remain 
invariant  with  respect  to  a  90-degree  material  rotation.  Similar  invariance  has  been 
previously  observed  by  Delale  and  Erdogan[18]  for  the  problem  of  an  orthotropic  strip 
with  internal  and  edge  cracks.  This  invariance  is  conditional  in  that  the  boundary 
conditions  at  x  =  0  axe  required  to  be  independent  of  the  y-coordinate[18].  In  the 
internal  and  edge  crack  problems  involving  membrane  stiffners,  where  the  boundary 
conditions  depend  on  the  y  coordinate,  it  is  shown  that  the  stress  intensity  factors 
indeed  vary  with  a  90-degree  material  rotation. 

From  its  derivation  it  is  cleax  that  the  integral  equation  (39)  is  valid  for  the  edge 
crack  as  well  as  for  the  intc  nal  crack  problem.  However,  for  c  =  0  it  is  seen  that  the 
kernel  N(s,  t)  becomes  unbounded  as  (x,  t)  — ►  0  simultaneously.  Thus,  apaxt  from  the 
Cauchy  singularity,  the  equation  contains  a  generalized  Cauchy  kernel  which  may  have 
an  effect  on  the  behavior  of  the  unknown  function  fx{t)  at  t  —  0.  This  behavior  may  be 
examined  quite  simply  by  following  a  function-theoretic  method  [19].  For  a  sectionally 
holomorphic  function  F(z)  defined  by  the  density  function  _/](£)  as 

F(z)  =  i  |  dt,  (46) 

o 

assuming  that 


A  (*)  =  -J  (7_  jp .  (0  <  Mb  -)  <  1), 

we  have  the  following  asymptotic  expression  [19] 


F(z)  = 


<M0) 


d?  sinj3~  ( z  —  cf  d ^  sinuiir  ( z  — 


+  F0(z), 


where  0\{t)  is  bounded  in  0  <  t  <  d  and  nonzero  at  t  =  0  and  t  =  d,  and  F0(z)  is 
either  bounded  or  has  singularities  that  are  of  lower  order  than  that  of  F(z).  If  we  now 
express  the  kernel  N(x,  t)  in  terms  of  the  following  partial  fractions 
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*>*«  -  cftrft  - 

_ 2  (2  +  P  Pa) 

^(px  ~  P2)  (*  -  P2z)  ’ 


2  (2  +  P  Px) 

CUPi  -  p2)  (*  -  ?!*) 


_  -p  ±  >|p2  -  4  _ 

Pi, 2  —  2  >  P  —  Q  +  2, 

the  integral  equation  (39)  becomes 


1  f  /,(*)  ,  ,  (p  +  2)  i  f  /,(() 


1  J  i  -  I 


(t  +  X) 


2  (2  +  P  V,)  i  f  /,(>) 
C*(P!  ~  P2)  ^  |  (<  -  Pxl) 


2  (2  +  P  P2)  1  4(t) 

QPi  -  P2)  77  (t  ~  P22) 


(0  <  X  <  d). 


Thus,  from  (48)  and  (51)  it  follows  that 


<*,(0)  J  ,  (?  +  2)  _  2(2  +  yy,)^T 

cFjP  [  Q  sin/3^-  (^(pj  —  p2)  szn/^Tr  p@ 

2(2  +  ^,)^  1  _  4,(4  _ 

c»(p,  -  Pj)  pfj  4?  (  4  -  xf~  C°  *  1  ' 

where  0(i)  represents  all  lower  order  terms.  Multiplying  (52)  by  ( d 
the  limit  as  x  — ►  d  we  obtain 


—  x)w  and  taking 


<M<0 


cotun  =  0. 


Similaxly,  multiplying  (52)  by  tP  and  letting  x-»0we  find 

<M0)  [  ,  (p  +  2)  _  2  (2  +  p  p, )  t#* 

^  \  C\  sin@K  Ch(pl  -  p2 )  sinfln  p@ 


2  (2  +  p  p2) 

ch{px  -  P2)  pf 


=  0. 
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Since  <^(0)  and  0x(tf)  are  non-zero,  from  (53)  and  (54)  we  obtain  the  following 
characteristic  equations  to  determine  u>  and  0. 


cotun  =  0, 


cot0ir  + 


( P  +  2) 
CA  sinfiir 


2  {2  +  p  p,}  * 

C»(P,  -  P,)  pf 


+ 


2  |2  +  p  p2j  e1^ 
C0pl  -  p2)  sinfin  p% 


0. 


(55) 


(56) 


From  (55)  it  is  seen  that  u  —  1/2  which  is  the  expected  result  for  the  singularity  at  a 
crack  tip  embedded  in  a  homogeneous  medium.  Equation  (56)  has  to  be  solved 
numerically  to  determine  the  value  of  0.  From  (50)  it  can  be  seen  that  px  and  p2  axe 
real  for  CA  >  0  and  are  complex  for  CA  <  0.  This  implies  that,  depending  on  the  sign  of 
CA,  (56)  gives  two  separate  characteristic  equations.  A  numerical  code  was  developed  to 
solve  (56)  for  positive  and  negative  values  of  C\  and  it  was  determined  that  it  does  not 
have  a  root  in  the  strip  0  <  Rt{0)  <1.  0  =  1  is  seen  to  be  a  root  but  is  unacceptable 
[19].  This  implies  that  the  unknown  function  jx(t)  does  not  have  a  power  singularity  at 
the  end  point  x  =  t  =  0. 

The  next  step  in  the  analysis  is  to  check  for  a  weaker,  namely  for  a  logarithmic 
singularity  in  the  unknown  function  when  c  =  0.  From  the  foregoing  analysis  we  note 
that  for  c  =  0  the  solution  of  the  integral  equation  is  of  the  form 


fi  (t)  =  !  ,  (57) 

(d-  t)i 

thus  from  (46)  and  (57)  the  asymptotic  behavior  of  F(z)  near  the  ends  z  =  0  and  z  —  d, 
respectively,  is  found  to  be 


F(z) 

F(z) 


0(0)  log  z 

— VW  +  F^}’ 
=JL  + 


(58) 

(59) 


where  d>{t),  Ft(z)  and  F2(z)  are  bounded  functions.  By  substituting  from  (57)  and  (58) 
into  (51)  it  may  be  easily  shown  that  the  coefficient  of  the  singular  term  log  x  vanishes. 
One  may,  therefore,  conclude  that  at  x  —  0  fi(x)  and,  consequently,  the  stresses  are 
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nonsingular. 


We  now  observe  that  for  c  =  0  by  defining  the  normalized  quantities 
i  =  i(r+i),  *  =  |(u+  1),  ^  =  <r(u),  (60) 

the  integral  equation  (39)  may  be  expressed  as 

1  /  V  i 

\  |  aT —  f  du  +  |  M3(r).  <(«)]  fx{t)  dt  =  —  1,  (  —  1  <  r  <  1),  (61) 

-l  -l 

where  N{s(r),  f(u)]  is  given  by  (40)  for  the  orthotropic  and  (42)  for  the  isotropic  half¬ 
plane.  From  (57)  it  is  clear  that  the  solution  of  (61)  is  of  the  form 

<r(u)  =  w(u)  G{u),  w(u)  =  (1  -  ttp(l  +  u)0.  (62) 


One  way  to  solve  the  problem  would  be  to  express  the  bounded  function  G(u )  in  terms 
of  a  series  of  the  related  orthogonal  polynomials,  in  this  case  the  Jacobi  polynomials,  as 
follows: 

00  I  1 

G(u)  =  (63) 

n= 0 

Equation  (61)  may  now  be  reduced  to  a  system  of  linear  algebraic  equations  in  An  by 
truncating  (63),  substituting  into  (61),  and  using  a  collocation  technique.  The 
collocation  points  are  chocen  according  to  (36).  The  following  expressions  are  useful  in 
evaluating  the  singular  integral[20]: 


i  jp^)(i-<)°(i  +  ^  *  =  C0(QI  (1  _  ,ni  +  3)v.«.  % 

-l 

Ja+p)  T(a)  T(n  -f  0  4-1)  ,  ... 

r(n  +  Q  +  /3  +  1)  F{(n+1)’  (  -  1  -  a  -  0),  (1  -  a),  (-J-5)}  , 


1  f  P(°-  0)(‘)  (1  -  1  f  (1  +  tf 
n  J  t  +  u 


(64) 

CO0 K  (1  +  U)a{l  -  U)V/’  °\u) 


(65) 

The  hypergeometric  function  F(n,  a,  b,  c),  appearing  in  (64)  and  (65),  is  computed  by 
using  the  summation  formula  given  in  [20].  The  integral  containing  the  generalized 


2(a+/J)  T{3)  T(n  4-  a  +1)  A/ 
T(n  +  a  +  0  +  l)  Fl( 
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Cauchy  kernel  is  evaluated  by  Gaussian  quadrature.  The  convergence  behavior  of  the 
numerical  scheme  is  presented  in  Table  3  which  gives  the  solution  for  the  isotropic  half¬ 
plane.  It  is  seen  from  the  table  that  excellent  convergence  is  obtained  with  an 
increasing  number  of  terms  in  the  series.  For  N  —  45  the  value  of  normalized  stress 
intensity  factor  matches  the  exact  value  computed  in  [21]  upto  seven  digits  beyond  the 
decimal  point.  The  exact  value  of  1.12152226  was  computed  in  [21]  where  it  was  stated 
that  the  only  possible  error  might  be  in  the  last  digit.  This  technique  is  used  in  all  edge 
crack  problems  in  this  study. 

After  determining  the  unknown  function,  the  stress  intensity  factor  for  an  edge 
crack  in  an  orthotropic  half-plane  may  be  obtained  from  (44).  For  different  Ch  values, 
the  normalized  stress  intensity  factor  is  tabulated  in  Table  4.  The  results  are  seen  to 
deviate  significantly  from  the  the  familiar  isotropic  result  (=  1.1215),  when  Ch  is 
different  from  zero.  As  in  the  internal  crack  case,  the  stress  intensity  factor  is  a 
decreasing  function  of  Ch.  Convergence  in  the  numerical  solution  was  observed  to  be 
very  rapid  for  positive  values  of  CK.  However,  it  becomes  increasingly  difficult  to  obtain 
convergent  results  for  CK  <  —  1.  The  results  presented  in  Table  4  axe  accurate  upto  the 
digits  shown.  They  have  been  obtained  by  increasing  the  number  of  terms  in  the  series 
and  the  quadrature  points.  Convergence  for  Ch  <  —  3.9  was  extremely  difficult  and  the 
present  scheme  yields  unreliable  results.  By  letting  CK  —  —  4  +  e,  a  singular 
perturbation  technique  may  be  developed  to  solve  the  problem.  However,  the  results 
would  be  of  no  practical  interest.  In  the  other  limiting  case,  for  Ch  — *•  oo  from  (27)  it 
may  be  shown  that  the  Fredholm  kernel  becomes 


N{x,  t) 


I  1 

z  (t  +  x)  ' 


(66) 


Thus,  by  substituting  from  (66)  into  (39)  and  assuming  /j(t)  =  —  the  integral 

equation  may  be  reduced  to 


i 


(-d<  x<  i), 


(67) 


giving 


From  (44)  and  (68)  it  then  follows  that 
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(69) 


*i(<i)  =  a.  "R , 

that  is,  in  Table  4  the  normalized  stress  intensity  factor  would  approach  one  as  CK  goes 
to  infinity. 

5.  Solution  of  the  Crack  Problem  for  a  Single  Stiffener 

Referring  to  Fig.l,  for  c  >  0,  the  system  of  singular  integral  equations  (15)-(17) 
may  be  solved  by  using,  for  example,  the  technique  described  in  the  previous  section  in 
a  straightforward  manner.  The  quantities  of  interest  to  be  calculated  are  the  stress 
intensity  factors  which,  for  the  case  of  a  single  film  shown  in  Fig.l  may  be  defined  as 
follows: 


*l(<0  = 

Urn  ^,j2(c  -  x)  crJS(x,  0)  , 

(70) 

*i(«0  = 

Urn  ^2(x  -  d)  cryy(x,  0)  , 

(71) 

*a(c)  = 

lim  J2(c  -  x)  <r  (x,  0)  , 

(72) 

fcj(d)  = 

lim  J2(x  -  d)  a^x,  0)  , 

(73) 

—  2a) 

=  y  2„  V  "b  ^a)  <Xty(y,  0)  , 

(74) 

*2(0)  = 

hmQ\-2y  (T^y,  0). 

(75) 

An  examination  of  the  kernels  N given  by  (18)-(26)  would  show  that,  for  a  fixed  crack 
length  (d  —  c)  and  the  film  length  2a,  as  the  crack  distance  (d  +  c)/2  goes  to  infinity, 
the  kernels  Nn,  Nl3,  N23,  N31  and  N32  would  go  to  zero  and  the  integral  equations  would 
uncouple.  Also,  since  jVn  and  N22  axe  invariant  with  respect  to  a  90-degree  material 
rotation,  the  crack  tip  stress  intensity  factors  would  remain  unchanged  through  such  a 
rotation.  However,  since  N33  is  explicitly  dependent  on  A  =  (Cu/C22)5,  even  in  the 
uncoupled  case  the  film  end  stress  intensity  factors  would  be  affected  by  a  90-degree 
material  rotation.  In  the  coupled  case  all  kernels  are  dependent  on  A.  Consequently, 
the  crack  tip  as  well  as  the  film  end  stress  intensity  factors  would  be  affected  by  a  90- 
degree  material  rotation. 

Some  sample  results  showing  the  normalized  stress  intensity  factors  in  an 
orthotropic  half-plane  stiffened  by  a  single  membrane  and  containing  an  internal  crack 
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axe  given  in  Tables  5-7.  In  these  examples,  too,  it  is  assumed  that  the  crack  angle  6 
shown  in  Fig.  1  is  zero  and  the  elastic  properties  of  the  substrate  are  given  in  Table  1. 
Note  that  the  problem  has  two  dimensionless  length  parameters,  {(d  —  c)/2a} 

characterizing  the  relative  crack  length  and  {(d  +  c)f 2a}  giving  the  relative  crack 
location.  Aside  from  the  magnitude  of  the  external  load  <70,  the  third  additional 
variable  in  the  problem  is  the  relative  compliance  of  the  film  characterized  by 


(76) 


Table  5  shows  the  results  for  a  fixed  crack  length  {(d  —  c)/2a}  =  1,  and  Table  6  for  a 
fixed  (average)  distance  {(d  4-  c)/2a}  =  1.  From  Table  5  it  may  be  seen  that  as  the 
stiffness  of  the  film  decreases,  (i.e.  as  X}  increases)  the  primary  crack  tip  stress  intensity 
factors  k\(c)  and  ^(d)  increase  and  approach  the  values  given  in  Table  2  (for  which  Xj  = 
oc).  As  the  crack  distance  from  the  surface  increases,  k^(c)  and  fc^d)  tend  to  zero  and 
the  remaining  stress  intensity  factors  approach  the  uncoupled  results  given  in  Table  2 
and  Fig. 3  regardless  of  the  value  of  the  film  stiffness  Xj.  The  rapid  increase  in  fc^O),  the 
stress  intensity  factor  for  the  interface  shear  at  the  film  end  x  =  0,  as  the  crack  tip 
approaches  the  surface,  (i.e.  as  c  —  0  or  {(d  +  c)/ 2a}  — ■  1)  is  due  to  the  interaction  of 
singular  stress  fields  at  the  crack  tip  and  the  film  end.  Table  6  shows  some  sample 
results  for  the  normalized  stress  intensity  factors  for  a  fixed  crack  distance  {(d  +  c)/2a} 
=  1  and  various  values  of  the  crack  length  {(d  —  c)/2a)  and  film  compliance  coefficient 
Xf.  Referring  to  the  integral  equation  (30)  and  Fig. 3,  note  that  the  normalizing  stress  in 
the  figure  is 


<?0 


A  <T0 

IQT!’ 


(77) 


which,  for  isotropic  materials  becomes  cr'0  —  <?0I2..  Therefore,  in  comparing  the  results 
in  Tables  5  and  6  with  that  of  Fig. 3,  the  values  given  in  Fig. 3  need  to  be  multiplied  by 
a  factor  of  (A/,JCa  +  4}. 

The  results  found  in  this  study  for  isotropic'  materials  compare  well  with  that  of 
[8]  (see  also  [17]).  The  slight  discrepancy  that  may  exist  is  due  to  different  techniques 
used  in  the  respective  numerical  analyses.  The  effect  of  90-degree  material  rotation  on 
the  stress  intensity  factors  is  shown  in  Table  7.  The  example  is  given  only  for  CA  = 
6.8222,  Material  1A.  Except  for  a  90-degree  rotation,  the  substrate  material  (Mat.lA) 
used  in  Table  7  is  the  same  as  that  used  in  Tables  5  and  6  (Mat.l),  and  in  both  cases  CA 
=  6.8222.  That  is,  referring  to  Table  1,  it  is  assumed  that  in  Tables  5  and  6  Cu  = 
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1.0100  x  10n,  C22  =  0.3592  x  1011  and  in  Table  7  Cn  =  0.3592  x  10u,  C22  =  1.0  1  00  x 
1011  N/rn2,  C12,  Cqq  and  CA  being  the  same.  Note  that  the  stress  intensity  factors  for 
the  material  orientation  used  in  Tables  5  and  6  axe  consistently  higher  than  that  given 
in  Table  7  and  the  differences  appear  to  be  significant. 

In  the  important  case  of  an  edge  crack  (c  =  0,  9  =  0,  Fig.l)  the  integral 
equations  (15)-(17)  are  still  valid.  However,  they  now  have  generalized  Cauchy  kernels 
which  become  unbounded  as  x,  t,  y,  y0  approach  the  end  point  x  =  0,  y  =  0 
simultaneously.  Unlike  the  edge  crack  problem  in  the  absence  of  a  stiffner,  described  in 
the  previous  section,  in  this  case  it  is  not  possible  to  make  an  intuitive  statement  about 
the  nature  of  the  singularity  at  the  common  end  points  of  the  crack  and  the  film.  This 
singularity  may  be  examined  in  a  straightforward  manner  by  defining 


T  AM  =  tr  ■  <°  <  Re M)  <  ■> 

.  (0  <  R<tf,  7)  <  1) 

A  t  (a  —  t)' 

fM  =  : — ■  (°  <  W,  c)  <  t) 
( -  y0y  (y0  +  2a)Q 


(78) 

(79) 

(80) 


and  by  following  a  function- theoretic  mathod  briefly  described  in  the  previous  section. 
For  the  remote  singular  points  the  characteristic  equations  are  found  to  be 


coLjk  =  0,  cot"fK  —  0,  cotan  =  0. 


(81) 


giving  a;  =  7  =  a  =  1/2.  The  characteristic  equation  found  for  the  endpoint  x  =  0  =  y 
is  rather  complicated[17].  However,  its  examination  shows  that  it  has  no  root  in  the 
complex  strip  0  <  Re(/3)  <  1.  Thus  the  unknown  functions  /x,  /2  and  /3  have  no  power 
singularity  at  x  =  y  =  0  and,  hence,  may  be  expressed  as 


Gi(t) 

TT-7 


c2(t) 


fziVo) 


GM 

Ty7+ 2a’ 


(82) 


where  G„  (t  =  1,  2,  3)  are  bounded  functions  and  are  nonzero  at  the  end  points. 
Searching  now  for  a  logarithmic  singularity  at  the  point  x  =  y  =  0,  after  performing  the 
appropriate  asymptotic  analysis  from  (15)-(17)  and  (82)  it  be  shown  that  the 
coefficients  of  the  logarithmic  terms  cancel  out  and 


-  C?3(0)  log  -  y  <  o c,  y  <  0, 


(83) 
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giving  G3(0)  =  0.  This  would  imply  that  the  apex  of  the  90-degree  wedge  formed  by 
the  intersection  of  the  crack  and  the  membrane  is  stress  free. 

In  the  case  of  a  single  stiffener  and  an  edge  crack  (Fig.  1 ,  c  =  0,  9  =  0),  after 
normalizing  the  intervals,  from  (82)  it  may  be  shown  that  for  fx  and  /2  the  weight 
function  and  the  related  orthogonal  polynomial  are  (1  —  t)  *  (1  +  t)°  and  °\t)  and 
for  /3  (1  —  t)°  (1  +  t)  2  and  *\t),  respectively.  The  technique  to  solve  the  integral 
equations  would  then  be  quite  similar  to  that  described  in  the  previous  section.  Note 
that  in  this  edge  crack  problem  the  equilibrium  condition  (13)  is  the  only  external 
condition  that  needs  to  be  satisfied  in  solving  the  integral  equations. 

After  solving  the  integral  equations  the  stress  intensity  factors  are  determined 
from  (71),  (73)  and  (74).  Some  sample  results  obtained  for  Ch  =  6.8222,  3.0358, 
—  1.0944  and  —  1.3137  corresponding  to  materials  1,  2,  5  and  6,  respectively  (see  Table 
1),  are  shown  in  Table  8  for  various  stiffness  ratios  \j  and  relative  dimensions  d/a. 
Physically  it  is  expected  that  as  (d/a)  — >  oo  (for  all  values  of  and  as  — »  oo  (  for  all 

values  of  (d/a)),  k? (d)  would  approach  zero  and  £j(d)  would  tend  to  the  edge  crack 
results  for  an  unstiffened  orthotropic  half-plane  given  in  Table  4.  This  trend  may 
indeed  be  observed  in  Table  8.  Similarly,  as  (d/a)  — *■  0  the  stress  intensity  factor  for 
the  interface  shear,  k^(  —  2a)  would  be  expected  to  approach  the  results  given  in  Fig.3. 
This,  too,  may  be  seen  in  Table  8  again  with  the  understanding  that  to  conform  to  the 
notation  given  in  Fig.3,  the  results  given  in  Table  8  must  be  multiplied  by 
J  (CK  +  4)/A.  Note  that  the  interaction  between  the  membrane  stiffener  and  the  crack 
is  much  stronger  in  the  edge  crack  case  than  it  is  for  the  internal  crack  problem.  The 
table  shows  that  generally  all  stress  intensity  factors  decrease  with  both  increasing  CA 
and  increasing  compliance  constant  Xj.  One  of  the  more  significant  results  of  the  edge 
crack  problem  is,  of  course,  the  fact  that  the  interface  shear  stress  intensity  factor  ^(0) 
at  the  film  end  x  =  0  is  zero.  The  effect  of  a  90-degree  material  rotation  on  the  stress 
intensity  factors  is  shown  in  Table  9.  The  results  given  in  Table  8  are  based  on  the 
material  constants  shown  in  Table  1.  In  Table  9  the  90-degree  rotation  is  accomplished 
simply  by  interchanging  Cn  and  C22  in  the  substrate.  In  this  case,  too,  the  stress 
intensity  factors  for  the  material  orientation  corresponding  to  Table  8  are  seen  to  be 
consistently  greater  than  that  given  in  Table  9.  Not  unexpectedly,  the  differences 
between  the  two  sets  of  results  are  seen  to  be  much  more  pronounced  for  small  values  of 
(d/a). 

Some  sample  results  for  the  distribution  of  the  interface  shear  stress  in  an 
isotropic  half-plane  containing  a  crack  and  stiffened  by  a  single  membrane  axe  also  given 
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in  Fig. 5.  Curve  1  shows  the  result  for  the  edge  crack  case  c  =  0.  Note  that  in  this 
particular  case  at  the  endpoint  y  =  0,  the  interface  shear  stress  is  zero.  For  reference, 
the  figure  also  shows  the  shear  stress  for  c  =  oo  (curve  4)  corresponding  to  the 
uncoupled  problem.  Also  note  that  as  the  ligament  size  c  decreases,  there  is  a 
considerable  increase  in  the  interface  shear  near  y  =  0. 

6.  The  Broken  Film  on  a  Cracked  Substrate 

The  problem  is  described  in  Fig. 2.  In  this  case  by  observing  that  /2(t)  =  0  and 
by  using  the  symmetry  conditions,  the  integral  equations  (15)-(17)  may  be  expressed  as 

d  ,,  ,  d  o 

i  j  dt  +  Qn(x>  *)  IM  dt  +  Q\s(x>  Vo)  h{Vo)  dVo 

0  0  -2a 

=  &(*),  (0  <  X  <  d),  (84) 

0  r  /  \  d  0 

if  Vo  -  y  dy°  +  &i( y ’  *)  /i(^)  dt  +  f  yo)  h(Vo)  dVo 

■  2a  0  -2a 

=  92(y),  ( ~2a  <  y  <  0)’  (85) 


where 


Qn(i,  t)  — 


QisU-  Vo)  = 


2  +  4) _ {2  x  -t-  t  {Ck  +  2)}  \ 

:r  Ch  \2  {x  +  t)  +  ( C7A  +  2)  x  ^  +  t2  }/  ’ 

4  jjg+ 1)  ^  A3 
Kn  {(y*  +  A2  2^)2  +  Ch  2^}  ’ 


Qn(y>  t)  =  - 


_ Kn(  Q  +  4)  t  A  V 

*  {(i2  +  A2  j/2)2  +  Ch  yW}  ’ 


Q23(y-  yo)  = 


i  _  P\2  *H(y  -  y0) 

X  (Vo  +  y)  P22  h](Ch+  4) 


(86) 

(87) 

(88) 

(89) 


In  this  case  too,  the  asymptotic  analysis  shows  that  at  the  end  point  x  =  y  =  0  the 
functions  /t  and  f3  are  bounded  and  have  the  form  as  given  by  (82)  (see  [17]).  The 
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solution  may,  therefore,  be  obtained  by  normalizing  the  intervals  and  by  expressing  Gx 
and  G2  in  terms  of  series  of  Jacobi  polynomials. 

Some  sample  results  for  normalized  stress  intensity  factors  are  listed  in  Table  10. 
Table  11  shows  the  effect  of  the  90-degree  material  rotation.  Compaxing  the  results 
given  in  Tables  8  and  9  with  those  of  Tables  10  and  11,  it  may  be  seen  that  the  stress 
intensity  factors  for  the  case  of  two  symmetric  membranes  are  consistently  greater  than 
the  single  membrane  case.  It  may  also  be  observed  that,  in  general,  the  stress  intensity 
factors  increase  with  decreasing  C 

7.  Conclusions 

The  membrane  model  used  in  this  paper  for  thin  film  overlays  on  elastic 
substrates  appears  to  be  a  good  first  approximation  to  various  stress  concentration 
problems  in  the  composite  medium.  The  results  obtained  in  the  form  of  stress  intensity 
factors  may  be  useful  in  the  study  of  such  mechanical  failure  problems  as  film 
debonding,  film  rupture  and  substrate  cracking.  The  analysis  was  carried  out  for 
orthotropic  films  and  orthotropic  substrates.  The  results  found  indicate  that  the 
material  orthotropv  may  have  a  significant  influence  on  the  stress  intensity  factors  at 
the  crack  tips  and  along  the  film  edges.  One  major  limitation  of  the  model  is  neglecting 
the  bending  stiffness  of  the  film.  As  a  consequence  of  this  the  normal  component  of  the 
stress  along  the  interface  is  assumed  to  be  zero,  which  is  not  realistic.  Since  the  elastic 
continuum  model  for  thin  films  leads  to  rather  insurmountable  analytical  and 
computational  difficulties,  a  “plate”  model  appears  to  be  the  obvious  compromise.  The 
crack  problem  for  an  orthotropic  half-plane  stiffened  by  a  plate  will  be  the  topic  of  our 
next  paper  on  the  subject. 
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Appendix  A 

Green’s  Functions  for  Stresses  in  an  Orthotropic  Half-Plane 


Consider  an  orthotropic  elastic  solid  having  the  constitutive  relations 


where  C1;  is  the  stiffness  matrix.  The  Green’s  functions  due  to  a  pair  of  dislocations 
at  a  point  (ij,  yx)  in  the  half-plane,  appearing  in  (1)  are  found  to  be  [12] 

K\\  {x,  y,  xv  y{)  = 

-  -  S’’  {(j  +  -  J2(s/-  y,)2}} 

+  M  (2  c  +  C‘ (l  +  ^  "  £(2(i  +  i-k 

+  C,{C,  i,  r  +  x  (x  +  z,)  (x  +  3  X,}  +  zr\2(y  -  Vi)2}}j .  (A2) 

Xu  U  y>  *n  Vi)  = 

-  {’Sr’ (<c*  +  3) (*  -  *>)2  + *2 (»  -  j.)2} 

-  ~ ~ -  —  {(C,  +  3)  (x  +  z,)2  +  A2  (y  -  Jli)2}j 

+  l2  (  +  c»  (*  +  *.)2} 
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{2  <  +  c,(a?  +  4}} , 


*21  (*>  y>  ®i.  yi)  = 

_  (s  ^  ft)  |(x  +  ij*  —  A2  (v  —  V,)1}! 

+  :r|{liLir!]{2<  +  c‘(I  +  1.)2} 

-  (2  c  +  ^  +  4}} . 

*22  y»  ^i7  Vi)  = 

ft  {— XT^  {<*  -  *■>*  -  A2  (S'  -  ft)2} 

-  (£ftal  {(*  +  *.)*  -  *2  (y  -  ft)2}} 

-  rev  f^r1  (2<  +  Ci  A‘ (y  "  *>’)  -  h  <2  <*  +  *)  < 

+  CK{Ch  X  ^  +  *1  (i  +  *,)  (3  I  +  *i)  +  *  A2(y  -  yi)2}}| , 

*3i  (x>  V’  *1.  yft  = 

ft  ((l  -  ^)2  +  (C*  +  3)  *2  (y  -  ft)2} 

~  " — {(x  +  3=1  (2  +  (^  +  3)  -S2  (y  -  ft)2}} 

-  #5V  {2  c  +  <7,  A*  (y  -  ft)-}  -  £  {2  (,  +  *,)  c 

+  Cj{C,  x  x?  +  x,  (x  +  x,)  (3  x  +  Xj)  +  i  A2( y  -  ft)2}}}  , 

*32  (*>  y,  ®i,  V\)  = 


2: 


-  {(*  +  ^  -  #  (v  -  s'.)2}} 

-  vth  fj4rl]  (2  <  +  c*  *(»  -  *>2}  -  Iirrl1  f2  < 

+  CA  {(C,  +  4)  a,  +  2  x,x  +  A2  (y  -  !/i)}}1  , 


where 


*  -  (&)! 


K  —  (T1  4-  r  \  f  ^66  ( 6*12  ~  Cn) 

‘  12  C‘2  |c2I  (Cl,  +  Cu  +  2  C66)  /  ’ 

Cl2  ~  {^11  ^22 P> 

K 

"  “  A2  ’ 

^  _  (C12  4-  C12)  (C12  ~  Qa  ~  2  Qe) 

A  "  Cia 

C  =  (*  +  2i)2  +  a2  (y  -  y\)\ 

Xi  =  {(1  -  ax)2  +  A2  (V  —  yt)2}  +  C,  (*  -  I,)2  A2  (y  - 
X2  =  <2  +  C»  (1  +  i,)2  A2  (y  -  y,)2, 


(A10) 


(All) 


( A12) 


(A13) 


(A14) 


(A15) 


xs  =  (C  +  Ci  i,1)2  +  Ci  (1  -  z,)2  a2  (y  -  »i)2- 


(Ate) 


Similarly  the  Green’s  functions  for  a  pair  of  concentrated  forces  acting  at  a  point  (0,  y0 ) 
on  the  boundary  axe  given  by  [13] 


Ki3{x,  y,  y0)  = 

_ A3  x*{y  -  y0 

ir{x*  +  \2(CK  +  2)3?(y 


~  y0 )2  +  A4  (y  -  y0)4}  ’ 


(A17) 
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Kl4(x,  y,  y„) 


_ Ax»j(Ct  +  4) _ 

it  {i*  +  A2  (CK  +  2)  i?  (y  -  y.)2  +  A*  (y  -  y„)4}  ’ 


K23(x,  y,  y0)  = 


A3x(y  -  y0)2  i(Ck  +  i) 


ir  {z4  +  A2  (CK  +  2)  x2  (y  -  yB)2  +  A4  (y  —  y0)4}  ’ 


K2a{x,  y,  y0)  = 

_ Ax2  (y  -  y0)l(Ch  +  t) _ 

7T  (z4  +  A2  ( Ch  +  2)  x2  {y  -  y0)2  +  A4  (y  -  y0)4}  ’ 


(A18) 


(A19) 


(A20) 


K33{x,  y,  y0)  = 


A3  (y  -  y0)3  m  +  T) 


*  {x4  +  A2  ( CK  +  2)  I2  (y  -  y0)2  +  A4  (y  -  y0)4}  ’ 


(A21) 


K34(x,  y,  y0)  = 


A  x  (y  —  y0)2J(C*  +  4) 


7T  {z4  +  A2  (Ck  +  2 )  X2  (y  -  y0)2  +  A4  (y  -  y0)4} 


(A22) 
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Fig.l.  Schematic  diagram  illustrating  the  loading  and  geometry  for  the 
problem  of  an  internally  cracked  half-plane  stiffened  by  a  film. 


Fig. 2.  Schematic  diagram  illustrating  the  loading  and  geometry  for  the 
problem  of  an  edge  cracked  half-plane  stiffened  by  symmetric  film. 
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Material 

Qi 

xlO11  4 

m 

Ci  2 

xlO10  4 
m 

Bra 

Q 6 

xlO10  4 

m 

ck 

1 

1.0100 

2.7430 

0.3592 

0.4905 

6.8222 

2 

0.5966 

0.6764 

0.1712 

0.5592 

3.0358 

3 

0.2097 

0.8989 

0.2097 

0.5993 

0.0000 

4 

0.2088 

0.9012 

0.2101 

0.5971 

-0.0022 

5 

1.6800 

6.6000 

1.6800 

8.4000 

-1.0944 

6 

1.1904 

5.3840 

1.1904 

5.9520 

-1.3137 

Table  1.  Stiffness  coefficients  Cl}  for  the  orthotropic  half-planes  used  in  the  examples. 
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as 

6.8222 

3.0358 

0.0000 

d  +  c 

*i(c) 

Kid) 

Kic ) 

Kid) 

fci(c) 

Kid) 

d  —  c 

K 

K 

K 

K 

K 

2.6753 

1.2436 

2.7540 

1.2646 

2.8820 

1.3009 

2.1411 

1.2127 

2.2048 

1.2321 

2.3081 

1.2661 

1.8925 

1.1922 

1.9485 

1.2106 

2.0394 

1.2430 

1.7402 

1.1768 

1.7910 

1.1943 

1.8740 

1.2253 

1.6348 

1.1643 

1.6816 

1.1810 

1.7584 

1.2108 

1.12 

1.5563 

1.1537 

1.5999 

1.1698 

1.6719 

1.1985 

1.14 

1.4952 

1.1447 

1.5360 

1.1602 

1.6040 

1.1879 

1.16 

1.4459 

1.1367 

1.4849 

1.1517 

1.5487 

1.1786 

1.18 

1.4052 

1.1298 

1.4415 

1.1441 

1.5027 

1.1703 

1.20 

1.3709 

1.12J4 

1.4053 

1.1373 

1.4637 

1.1626 

1.24 

1.3164 

1.1126 

1.3473 

1.1254 

1.4008 

1.1493 

1.28 

1.2794 

1.1034 

1.3028 

1.1155 

1.3521 

1.1382 

1.32 

1.2418 

1.0954 

1.2677 

1.1069 

1.3131 

1.1284 

1.36 

1.2152 

1.0887 

1.2387 

1.0994 

1.2811 

1.1199 

1.1932 

1.0827 

1.2148 

1.0929 

1.2544 

1.1123 

1.1520 

1.0706 

1.1799 

1.0793 

1.2035 

1.0967 

1.1236 

1.0611 

1.1387 

1.0690 

1.1676 

1.0844 

1.0872 

1.0475 

1.0984 

1.0537 

1.1203 

1.0664 

1.0653 

1.0382 

1.0740 

1.0434 

1.0913 

1.0539 

— +oo 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

Td.ole  2.  Normalized  stress  intensity  factors  in  an  orthotropic  half-plane  with  an 
internal  crack,  Fig. 4,  9  =  0,  k0  =  a oy  (d  —  c)/2,  Materials  1-3. 
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k(d) 

7J7 

N 

1.121518230454 

10 

1.121522319145 

15 

1.121522334904 

20 

1.121522287226 

25 

1.121522267425 

30 

1.121522259954 

35 

1.121522255943 

45 

Table  3.  Convergence  of  the  calculated  stress  intensity  factor  for  an  edge  crack  in  an 
isotropic  half-plane  under  tension. 
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d  +  c 


2a 


M<0 


d- 


2 


f  =  0.00 


1.1 

1.0778 

1.4322 

0.0454 

-0.0965 

0.9847 

-0.4760 

1.5 

1.0288 

1.0416 

0.0309 

0.0184 

0.5718 

-0.5219 

2.0 

1.0170 

1.0132 

0.0223 

0.0269 

0.4758 

-0.4946 

3.0 

1.0089 

1.0073 

0.0136 

0.0188 

0.4222 

-0.4466 

5.0 

1.0040 

1.0038 

0.0065 

0.0088 

0.4006 

-0.4123 

A;  =  0.20 

1.1 

1.0823 

1.4367 

0.0438 

-0.0882 

0.9495 

-0.4488 

1.5 

1.0313 

1.0481 

0.0291 

0.0178 

0.5440 

-0.4955 

2.0 

1.0133 

1.0165 

0.0209 

0.0253 

0.4521 

-0.4708 

3.0 

1.0094 

1.0084 

0.0127 

0.0176 

0.4016 

-0.4256 

5.0 

1.0042 

1.0041 

0.0061 

0.0083 

0.3816 

-0.3930 

=  1.00 

1.1 

1.0961 

1.4521 

0.0386 

-0.0634 

0.8382 

-0.3678 

1.5 

1.0390 

1.0679 

0.0238 

0.0158 

0.4592 

-0.4157 

2.0 

1.0223 

1.0263 

0.0168 

0.0205 

0.3803 

-0.3989 

3.0 

1.0110 

1.0116 

0.0101 

0.0141 

0.3394 

-0.3621 

5.0 

1.0046 

1.0048 

0.0048 

0.0066 

0.3239 

-0.3346 

Table  5.  Stress  intensity  factors  in  an  orthotropic  half-plane  stiffened  by  a  single 
membrane  and  containing  an  internal  crack;  Fig.l,  B  —  0,  {(d  —  c)/2o}  =  1,  = 

{l3lJi322){a/h),  Material  1,  Ch  =  6.8222. 
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k,(d) 

M<0 

Md) 

^(c) 

f  =  0.00 


0.10 

0.9465 

0.9401 

0.0193 

0.0172 

0.3962 

-0.3952 

0.25 

0.9560 

0.9417 

0.0210 

0.0151 

0.4099 

-0.4033 

0.50 

0.9806 

0.9677 

0.0253 

0.0070 

0.4684 

-0.4293 

0.75 

1.0212 

1.0859 

0.0336 

0.0235 

0.6257 

-0.4582 

0.9° 

1.0685 

1.3993 

0.0431 

-0.0997 

0.9256 

-0.4536 

A,  =  0.20 


0.10 

0.9500 

0.9440 

0.0182 

0.0163 

0.3775 

-0.3766 

0.25 

0.9594 

0.9461 

0.0198 

0.0144 

0.3905 

-0.3841 

0.50 

0.9842 

0.9736 

0.0240 

0.0071 

0.4463 

-0.4081 

0.75 

1.0254 

1.0930 

0.0320 

-0.0207 

0.5983 

-0.4339 

0.90 

1.0731 

1.4033 

0.0417 

-0.0913 

0.8923 

-0.4280 

.9556 

.9594 

.9912 

.1151 


,f  =  1.00 


U.0149 

0.0162 

0.0198 

0.0273 

0.0371 


0.0135 

0.3210 

-0.3202 

0.0122 

0.3318 

-0.3260 

0.0071 

0.3791 

-0.3442 

-0.0128 

0.5134 

-0.3611 

-0.0663 

0.7874 

-0.3514 

Table  6.  Stress  intensity  factors  in  an  orthotropic  half-plane  stiffened  by  a  single 
membrane  and  containing  an  interned  crack;  Fig.l,  6  =  0,  {( d  +  c)/ 2a}  =  1,  Ay  = 
(. 3l22/022)(a/h),  Material  1,  CK  =  6.8222. 
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d  +  c 


2a 


kiid) 


d—  c 


2 


\{  =  0.00 


M<0  I  MO) 


d- 

2 


1.3985 

0.0276 

-0.0174 

0.6571 

-0.4018 

1.0817 

0.0139 

0.0196 

0.3452 

-0.3744 

1.0424 

0.0080 

0.0136 

0.2789 

-0.3161 

1.0186 

0.0039 

0.0065 

0.2476 

-0.2678 

1.0068 

0.0014 

0.0022 

0.2372 

-0.2435 

..4043 

0.0268 

-0.0159 

0.6409 

-0.3880 

..0843 

0.0134 

0.0189 

C.3343 

-0.3634 

L.0433 

0.0077 

0.0130 

0.2702 

-0.3071 

1.0188 

0.0037 

0.0062 

0.2402 

-0.2602 

1.0068 

0.0014 

0.0021 

0.2304 

-0.2365 

A  f  =  1.00 


1.4244 

0.0240 

-0.0110 

0.5852 

-0.3420 

L  .0932 

0.0116 

0.0166 

0.2974 

-0.3261 

1.0463 

0.0067 

0.0113 

0.2411 

-0.2767 

1.0195 

0.0032 

0.0054 

0.2155 

-0.2347 

1.0070 

0.0012 

0.0018 

0.2073 

-0.2132 

Table  7.  The  effect  of  a  90-degree  material  rotation  on  the  stress  intensity  factors  in  an 
orthotropic  half- plane  stiffened  by  a  single  membrane  and  containing  an  internal  crack; 
Fig. 5,  0  =  0,  {{d  -  c)/2a}  =  1,  \f=  {d\j022)(a/h^  Material  1A,  Ch  =  6.8222. 
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h 

d 

a 

<?o  fd 

k?{d) 

<?0  Vd 

cK 

=  6.8222 

0.10 

1.6923 

-0.2201 

0.0 

0.25 

1.2998 

-0.0973 

0.50 

1.1138 

-0.0356 

1.00 

-0.0005 

0.10 

1.6568 

-0.2078 

-0.3386 

0.2 

0.25 

1.2492 

-0.0916 

-0.3051 

0.50 

1.1114 

-0.0336 

-0.2333 

1.00 

1.0702 

-0.0005 

-0.1221 

0.10 

1.5480 

-0.1701 

-0.2923 

1.0 

0.25 

1.2159 

-0.0742 

-0.2677 

0.50 

1.1035 

-0.0272 

-0.2081 

i 

1.00 

1.0695 

-0.0004 

-0.1115 

Table  8.  Stress  intensity  factors  in  an  ortt 
membrane  and  containing  an  edge  crack;  Fi 
Materials  1  and  2. 
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*i( d ) 

-2a) 

a0  -fd 

<7,  Vd 

■H 

Ck 

=  3.0358 

-0.4747 

-0.4345 

-0.0671 

-0.3515 

-0.0134 

-0.2059 

1.8592 

-0.3118 

-0.4485 

1.3593 

-0.1482 

-0.4130 

1.1690 

-0.0622 

-0.3364 

1.0975 

-0.0124 

-0.1989 

1.6869 

-0.2439 

-0.3726 

1.2985 

-0.1143 

-0.3492 

1.1514 

-0.0478 

-0.2903 

1.0954 

-0.0100 

-0.1767 

jtropic  half-plane  stiffened  by  a  single 

.1,  6  =  0,  c  =  0,  A,  =  (P\2/Pn)(a/h), 


1 

d 

a 

M<9 

fcj(d) 

TJd 

—  2a) 

*,(<0 

~a~Td 

e 

CM 

1 

<y0  ^5 

KH 

cr0  {d 

m 

=  - 1.0944 

Ck  =  -1.3137 

0.10 

2.1371 

-0.3809 

-0.5481 

2.1853 

-0.4030 

-0.5708 

0.0 

0.25 

1.5102 

-0.1818 

-0.4938 

1.5377 

-0.1940 

-0.5159 

0.50 

1.2592 

-0.0742 

-0.3941 

1.2751 

-0.0805 

-0.4145 

1.00 

1.1585 

-0.0110 

-0.2030 

1.1668 

-0.0124 

-0.2168 

2.0570 

-0.3508 

-0.5140 

2.0990 

-0.3701 

-0.5340 

m 

1.4814 

-0.1670 

-0.4666 

1.5061 

-0.1776 

-0.4864 

1.2509 

-0.0684 

-0.3760 

1.2659 

-0.0739 

-0.3946 

1.1578 

-0.0102 

-0.1964 

1.1659 

-0.0116 

-0.2095 

1.8369 

-0.2678 

-0.4190 

1.8645 

-0.2804 

-0.4325 

1.0 

1.4006 

-0.1256 

-0.3883 

1.4187 

-0.1325 

-0.4023 

1.2270 

-0.0517 

-0.3216 

1.2393 

-0.0554 

-0.3357 

1.1556 

-0.0008 

-0.1757 

1.1635 

-0.0009 

-0.1866 

Table  8.  Contd.,  Materials  5  and  6. 


\ 

d 

k(d) 

Jfcj(-2  a) 

fc 2(  —  2a) 

xi 

a 

<?„  ^ 

cr0  4d 

cr0 

°0  Jd 

m 

Ck 

=  6.8222 

cK 

=  3.0358 

0.10 

1.4105 

-0.0865 

-0.2042 

1.5155 

-0.1121 

-0.2443 

0.0 

0.25 

1.1389 

-0.0289 

-0.1578 

1.1889 

-0.0400 

-0.1941 

0.50 

1.0754 

-0.0006 

-0.0937 

1.1017 

-0.0008 

-0.1194 

1.00 

1.0658 

0.0000 

-0.0219 

1.0864 

0.0000 

-0.0264 

0.10 

1.3987 

-0.0836 

-0.1993 

1.4983 

-0.1076 

-0.2374 

m 

0.25 

1.1365 

-0.0279 

-0.1546 

1.1850 

-0.0384 

-0.1894 

0.50 

1.0752 

-0.0006 

-0.0921 

1.1012 

-0.0008 

-0.1171 

1.00 

1.0658 

0.0000 

-0.0218 

1.0865 

0.0000 

-0.0262 

0.10 

1.3586 

-0.0733 

-0.1824 

1.4415 

-0.0926 

-0.2140 

1.0 

0.25 

1.1284 

-0.0244 

-0.1433 

1.1721 

-0.0330 

-0.1735 

0.50 

1.0742 

-0.0005 

-0.0865 

1.0995 

-0.0007 

-0.1091 

1.00 

1.0658 

0.0000 

-0.0213 

1.0865 

0.0000 

-0.0257 

Table  9.  The  effect  of  a  90-degree  material  rotation  on  the  stress  intensity  factors  in  an 
orthotropic  half-plane  stiffened  by  a  single  membrane  and  containing  an  edge 
crack;Fig.l,  9  =  0,  c  =  0,  \}  =  (/3l22/ 022)(a/h),  Materials  1A  and  2A. 
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d 

a 

Kid) 

<y0  fd 

m 

Kid) 

<y0 

Kid) 

m 

1 

0.10 

2.4567 

-0.3718 

2.9870 

-0.5085 

2.9926 

-0.4978 

0.0 

0.25 

1.4614 

-0.3236 

1.6926 

-0.4475 

1.7074 

-0.4218 

0.50 

1.1613 

-0.2413 

1.2629 

-0.3531 

1.2809 

-0.3195 

1.00 

1.0750 

-0.1247 

1.1099 

-0.2051 

1.1360 

-0.1535 

0.10 

2.3702 

-0.3545 

2.8327 

-0.4771 

2.8477 

-0.4686 

0.2 

0.25 

1.4392 

-0.3112 

1.6480 

-0.4244 

1.6666 

-0.4014 

0.50 

1.1564 

-0.2335 

1.2510 

-0.3379 

1.2710 

-0.3070 

1.00 

1.0746 

-0.1216 

1.1085 

-0.1982 

1.1352 

-0.1494 

0.10 

2.1112 

-0.3021 

2.4076 

-0.3891 

2.4418 

-0.3852 

1.0 

0.25 

1.3705 

-0.2720 

1.5200 

-0.3564 

1.5475 

-0.3405 

0.50 

1.1409 

-0.2083 

1.2160 

-0.2914 

1.2410 

-0.2683 

1.00 

1.0732 

-0.1U0 

1.1043 

-0.1763 

1.1330 

-0.1357 

Table  10.  Stress  intensity  factors  in  an  orthotropic  half-plane  stiffened  by  two 
symmetric  membranes  and  containing  an  edge  crack;  Fig.2,  A j  =  fi22)(a/h), 

Materials  1-3. 
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h 

d 

a 

cr0  id 

!%( -  2a) 

k,(d) 
cr0  'Q 

■ 

K  (<*) 

<y0  ^ 

m 

ck  =  - 

-  1.0944 

BSE 

0.10 

2.9899 

-0.4882 

3.3935 

-0.5885 

3.5018 

-0.6139 

0.0 

0.25 

1.7063 

-0.4174 

1.8914 

-0.5077 

1.9418 

-0.5310 

0.50 

1.2805 

-0.3157 

1.3697 

-0.3951 

1.3947 

-0.4157 

1.00 

1.1359 

-0.1522 

1.1705 

-0.2022 

1.1804 

-0.2159 

0.10 

2.8453 

-0.4596 

3.1908 

-0.5477 

3.2821 

-0.5697 

0.2 

0.25 

1.6656 

-0.3975 

1.8310 

-0.4787 

1.8753 

-0.4994 

0.50 

1.2705 

-0.3037 

1.3534 

-0.3769 

1.3765 

-0.3958 

1.00 

1.1352 

-0.1482 

1.1691 

-0.1957 

1.1787 

-0.2088 

0.10 

2.4405 

-0.3780 

2.6537 

-0.4378 

2.7079 

-0.4520 

1.0 

0.25 

1.5469 

-0.3378 

1.6628 

-0.3958 

1.6930 

-0.4102 

0.50 

1.2407 

-0.2661 

1.3066 

-0.3224 

1.3244 

-0.3367 

1.00 

1.1329 

-0.1348 

1.1649 

-0.1752 

1.1738 

-0.1861 

Table  10.  Contd.,  Materials  4-6. 
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d 

a 

0.10 

0.0 

0.25 

0.50 

1.00 

0.10 
1.0  0.25 


fct( d )  fc;(  -2a)  Ud)  &;(  -2a)  fet( d !)  ^(-20) 


1.7900 

1.2115 

1.0850 

1.0658 


1.7639 

1.2068 

1.0844 

1.0658 


1.6761 

1.1907 

1.0825 

1.0659 


<7  0  ^ 


6.8222 

wm 

-0.2123 

1.9944 

-0.2551 

2.9962 

-0.4988 

-0.1585 

1.2909 

-0.1954 

1.7090 

-0.4227 

-0.0932 

1.1169 

-0.1190 

1.2816 

-0.3204 

-0.0216 

1.0865 

-0.0264 

1.1361 

-0.1542 

-0.2069 

1.9561 

-0.2473 

2.8507 

-0.4694 

-0.1552 

1.2832 

-0.1906 

1.6680 

-0.4022 

-0.0917 

1.1158 

-0.1167 

1.2716 

-0.3078 

-0.0215 

1.0865 

-C.0262 

1.1354 

-0.1501 

-0.1884 

1.8304 

-0.2215 

2.4437 

-0.3858 

-0.1439 

1.2574 

-0.1745 

1.5484 

-0.3410 

-0.0862 

1.1124 

-0.1088 

1.2414 

-0.2689 

-0.0210 

1.0866 

-0.0256 

1.1331 

-0.1363 

Table  11.  The  effect  of  a  90-degree  material  rotation  on  the  stress  intensity  factors  in 
an  orthotropic  half-plane  stiffened  by  two  symmetric  membranes  and  containing  am 
edge  crack;  Fig.2,  A,=  ({3l22/ 022)(a/h),  Materials  1A-3A. 
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LIST  OF  FIGURES 


Fig.l.  Schematic  diagram  illustrating  the  loading  and  geometry  for  the 
problem  of  an  internally  cracked  half-plane  stiffened  by  a  film. 

Fig. 2.  Schematic  diagram  illustrating  the  loading  and  geometry  for  the 

problem  of  an  edge  cracked  half-plane  stiffened  by  symmetric  films. 

Fig. 3.  Variation  of  the  normalized  film  edge  stress  intensity  factor 
with  stiffness  parameter  \t. 

Fig. 4.  Typical  variation  of  the  normalized  film  stress  along  the  contact 
region. 

Fig. 5.  Interface  shear  distribution  for  an  isotropic  half-plane  containing  a  crack  and 
stiffened  by  a  single  membrane  (Fig.  1,0  =  0);  the  curves  1  through  4 
correspond  to  c/a  =  0  (the  edge  crack),  0.1,  1.0  and  oo. 
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The  effect  of  a  90-degree  rotation  on  the  stress  intensity  factors  in  an 
orthotropic  half-plane  stiffened  by  a  single  membrane  and  containing 
an  internal  crack;  Fig.l,  5  =  0,  {(d  —  c)/ 2a}  =  1, 

A /=  (/4//y(aA),  Material  1A,  CK  =  6.8222. 

Stress  intensity  factors  in  an  orthotropic  half-plane  stiffened  by  a 
single  membrane  and  containing  an  edge  crack;  Fig.l,  5  =  0,  c  =  0, 

A f=  (0l22/022)(a/h),  Materials  1,  2,  5  and  6. 

The  effect  of  a  90- degree  rotation  on  the-  stress  intensity  factors  in  an 
orthotropic  half-plane  stiffened  by  a  single  membrane  and  containing 
an  edge  crack;  Fig.l,  5  =  0,  c  =  0,  A^  =  {022/ 022)(a/h),  Materials 
1A  and  2A. 

Stress  intensity  factors  in  an  orthotropic  half-plane  stiffened  by  two 
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symmetric  membranes  and  containing  an  edge  crack;  Fig. 2,  \j  = 
(P\2/022)(a/h),  Materials  1-6. 

Table  11.  The  effect  of  a  90-d^gree  material  rotation  on  the  stress  intensity 
factors  in  an  orthotropic  half-plane  stiffened  by  two  symmetric 
membranes  and  containing  an  edge  crack;  Fig. 2,  A j  = 
fl22)(a/h),  Materials  1A-3A. 
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